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Abstract 

We report measurements of the in-plane magnetic penetration depth AA(T) in single crystals 
of ErNi2B2C down to ~0.1 K using a tunnel-diode based, self-inductive technique at 21 MHz. 
We observe four features: (1) a slight dip in AA(T) at the Neel temperature T/v = 6.0 K, (2) a 
peak at T\yfm = 2.3 K, where a weak ferromagnetic component sets in, (3) another maximum at 
0.45 K, and (4) a final broad drop down to 0.1 K. Converting to superfluid density p s , we see that 
the antiferromagnetic order at 6 K only slightly depresses superconductivity. We seek to explain 
some of the above features in the context of antiferromagnetic superconductors, where competition 
between the antiferromagnetic molecular field and spin fluctuation scattering determines increased 
or decreased pairbreaking. Superfluid density data show only a slight decrease in pair density in 
the vicinity of the 2.3 K feature, thus supporting other evidences against bulk ferromagnetism in 
this temperature range. 



The magnetic members of the rare-earth (RE) nickel borocarbide family, REM2B2C 
(RE = Ho, Er, Dy, etc) have attracted much interest due to the interplay between mag- 
netism and superconductivity. ErNi 2 B 2 C, in particular, is a good candidate for study: 
superconductivity starts at T c « 11 K, before antiferromagnetic (AF) order sets ir*i at T/v 
m 6 K. In the AF state the Er spins are directed along the 6-axis, forming a transversely 
polarized, incommensurate spin- density- wave (SDW) state, with modulation vector modu- 
lation vector— 5 = 0.553a* (a* = 2n/a), before squaring up at lower temperatures 3 . Below 
7Vfm=2.3 K a net magnetization appears, superposed on a modulation with a periodicity 
of 20a*, confirming the microscopic existence of spontaneous weak ferromagnetism (WFM) 
with superconductivity 3 -^. 

The fact that T/v < T c enables us to study the influence of magnetism on superconductiv- 
ity. In particular, in this paper we study, via the penetration depth, the pairbreaking effects 
of the various magnetic orders on the superfluid response of this material. There have been 
several previous penetration depth measurements on ErNi2B2C. Jacobs et a/.— measured the 
microwave surface impedance of single- crystal ErNi2B2C from T c down to 4 K, but did not 
see the AF transition at 6 K. They concluded that the AF transition is not accompanied by 
changes in pairbreaking in zero field. Andreone et al. measured the microwave properties of 
ErNi2B2C thin films — microwave surface resistance 6 down to 2 K, and the change in pen- 
etration depth 7 from 2-5 K. They too, saw no feature at T/v, and attributed its absence to 
the smearing of the susceptibility x{^)- Gammel et al. performed small-angle neutron mea- 
surements (SANS) 8 on single-crystal ErNi2B2C down to 4 K, finding a decrease of A below 
Tn that they could not account for quantitatively. In this paper we present high-precision 
measurements of the in-plane magnetic penetration depth of single-crystal ErNi 2 B 2 C down 
to 0.1 K. We see features at T/v and Twfm, and ascribe these to the pairbreaking effects 
of AF order at T N = 6 K and the weak ferromagnetic ordering at T WFM = 2.3 K. We also 
observe a peak at 0.45 K, which we attribute to the presence of a spontaneous vortex phase 9 
(SVP), expected to occur in superconductors where ferromagnetism and superconductivity 
coexis t 10 ^ 11112 . The superfluid density graph indicates that these three magnetic orderings 
coexist with superconductivity, i.e. they do not destroy superconductivity in this material. 

Various theories of antiferromagnetic superconductors have been propose d 13 i 14 i 15 i 16 i 17 i 18 i 19 . 
We shall follow that of Chi and Nagi 1,9 , which is an extension of the mean-field model by Nass 
et a/ . 17 ' 18 to the regime where the superconducting gap A is finite, and it includes the effects 
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of spin-fluctuations, molecular field and impurities. In the Chi-Nagi (CN) model, which ap- 
plies specifically to superconductors with Tn < T c , two temperature regimes are separately 
considered. First, in the paramagnetic regime (T N < T < T c ), the depression of T c with 
respect to the non-magnetic counterparts, LuNi 2 B 2 C or YNi 2 B 2 C, is due to the exchange 
scattering of the conduction electrons from the spins of the RE Er ions. Assuming that the 
exchange interaction is weak, this paramagnetic phase of ErNi 2 B 2 C can be accounted for by 
the Abrikosov-Gorkov (AG) pairbreaking theory 29 . Second, in the AF phase (T < Tn), the 
effect of pairbreaking depends on the competition^ between the temperature-dependent AF 
molecular field (with parameter Hq(T)) and spin- fluctuation scattering of the conduction 
electrons, the latter by both magnetic RE ions (parameter l/r^ ) and non- magnetic impu- 
rities (parameter 1/Vi). The molecular field opens AF gaps on parts of the Fermi surface 
(FS), hence destroying the superconducting gap in those areas. The non-magnetic impurities 
do not affect the BCS state for an s-wave superconductor—, but weaken the effect of the AF 
field by destroying the pairing state for charge density waves or spin density waves 2 ^ — thus 
non-magnetic impurities promote the recovery of superconductivity. Moreover, the effect of 
the molecular field and spin fluctuations is governed by a sum rule^i, and the competition 
between them determines whether the AF phase gives increased or decreased pairbreaking 
below T/v- The total electronic effective magnetic scattering rate I/t^ is temperature- 
dependent and decreases with decreasing temperature (as the magnetic moments become 
more and more frozen). The assumptions of the CN model are: (1) the effect of inelastic 
scattering, which is relevant only for T <C T/v, can be ignored; (2) BCS s-wave pairing and 
(3) a one-dimensional (1-D) electron band that satisfies the nesting condition = — €k+Q. 

The following equations of the CN model were used^. The temperature-dependence of 
the superconducting gap is determined from 

(AG equation) .„ (|) = * (I) - 4, (I + , (1) 

(Renormalized frequency) u n ± = u n + Kp— + Kfc—^— , (2) 

2A+ 2A_ 

(Renormalized gap) A n ± = A ± H Q (T) + X T ^± + X ± ^p, (3) 
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(Gap equation) 
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where T c o is the superconducting transition temperature of the non-magnetic member of 
the borocarbide family LuN^P^C or YN^E^C, ip is the digamma function, X± and Y± are 
linear combinations of the magnetic (l/rf ), non-magnetic (1/n) and spin-orbit (l/r so ), 
scattering rates 
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A± = Pn± + A 2 ^] 1 / 2 ; ^n± = ^Vt/A n ±; = 7rT(2n+l) is the Matsubara frequency. In the 
paramagnetic phase the distinction between + and — is lost, giving uj n ± = ui n , A n ± = A n 
and so U n = u n /A n . 

The temperature-dependence of the superfluid density p s is given by^ 
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^^K^ + A^^^H^L+AL) 1 / 2 !. (9) 
In the paramagnetic phase (T N < T < T c ), p s is given by (P: paramagnetic) 



2 +A 2 a/2 



n>0 V 



(10) 



where £ 



-A 2 n ) 1/2 |. Note that this expression for the superfluid density is for materials 



in the dirty limit. This is consistent with the Er ions in ErNi 2 B 2 C being the "impurity" ion 
when compared to either LUM2B2C or YN^P^C. 



4 



We turn next to the parameters of the model. The effective magnetic scattering rate 
(I/T2 ) from RE ions (1/V.f) and magnetic impurities (1/t|) is given by 





~2 ~2 



(T > T N ) 
(T)) (T<T N ) 



(11) 



where — 

T 2 



27in R N(0)J(J+l)(gj-l) 2 I 2 . 



(12) 



The AF molecular field is given by 



H Q (T) = Hq(0)F(T) 



(13) 



where H Q (0) = n R I \ Qj - 1 | y/j(J+l). 



(14) 



np> is the concentration of RE ions, I is the exchange energy, gj is the Lande factor, and J 
is the total angular momentum of the RE ion. The function F(T) can be approximated by 
the empirical relation 



where v is a parameter obtained by fitting F(T) to sublattice magnetization data. 

The values of the renormalized frequencies u n ± and gaps A n± are determined self- 
consistently: for a fixed temperature T and Matsubara index n, one determines u n ± and 
A n ± from Equations 12 and El such that they also satisfy Eqn. |U After computing the w„±'s 
and A n ±'s for a fixed T, one then substitutes these values into Eqn. 0or E3to obtain the 
superfluid density p s at that temperature T. 

Details of sample growth and characterization are described in Ref. Q. The samples 
were then annealed according to conditions described in Ref. |3 The superconducting 
transitions of our sample were measured by low-field (if =5 G) magnetization, zero-field 
resistivity and zero-field specific-heat measurements. From magnetization data, the onset 
of superconducting diamagnetism appears at T=11.0 K and 90% of the full diamagnetic 
magnetization is reached at T=9.6 K. Resistivity data show a superconducting onset at a 
higher temperature of 11.3 K and zero resistivity at 9.6 K. The mid-point of the specific- 
heat jumpSi yields a T c of 10.1 K. A comparison of the three measurements show that bulk 
superconductivity occurs at T c w 10 K, whereas the initial decrease of resistivity at ~11 K 
may be due to some sort of filamentary superconductivity. 




(15) 
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The parameters of this model are determined as follows. We denote Ao and A(0) to be the 
zero-temperature superconducting gap amplitude of YNi 2 B 2 C and ErNi 2 B 2 C, respectively. 
Tunneling measurements^ yield A = 1.83T c0 . From the experimental values of T c (10.1 K 
for ErNi 2 B 2 C) andT c0 (15.5 K for YNi 2 B 2 C), Eqn. Ogives l/r 2 e// A = 0.227. We assume 1/t* 
= 0, in which case Eqn. 1111 gives l/r 2 R A =0.227. For ErNi 2 B 2 C, using the values n#=l/6, 
J=7.5, g-j=1.2, iV(0)=0.36 states/eV- atom-spin 27 , we obtain 7=0.024 eV from Ean.lT2lwhich 
is comparable with the experimental value of 0.031 eV— . This justifies our assuming l/r 2 
= 0, as any finite t% would make I even smaller than the experimental value. Equation IT41 
then gives Hq(0)/Ao=2.6. From the temperature-dependence of the magnetic Bragg peak 
intensity*^ below T/v we obtain v = 4.8 in Eqn. ED The only remaining free parameter 
of the theory is 1/r, the non- magnetic scattering rate, defined to be 1/r = 1/tl + 2/3r so . 
In an AFSC, usually 1/r ^> l/r^ 1 -, so for the present study we assume l/r so =0, such that 

1/T =1/71. 

To see the pairbreaking effects of the various magnetic orders we need to convert AA(T) 
to p s (T), the superfluid density. To determine p s (T) we need the value of A(0), which has 
been reported over a ranged from 700 A to 1150 A. We take A(0) to be a parameter in our 
model, keeping in mind that it has to be in the vicinity of the above two values. 

Figure H] shows the temperature-dependence of the in-plane penetration depth AA(T). 
The onset of superconductivity, T*, is 11.3 K, showing that this is a high-quality single 
crystal. We also see the following features: (1) a slight dip in AA(T) at T/v=6.0 K, (2) a 
peak at 7Vfm=2.3 K, (3) another maximum at 0.45 K, and (4) a final broad drop down 
to 0.1 K. We attribute the last two features to the presence of the SVP. For dirty AF 
superconductors, the penetration depth is expected to decrease below T N by both^ the 
susceptibility (x) and mean free path (/) as A ~ X' L / a/1 + 477%, where X' L ~ A^(l + £o/l). 
Neither effect, however, explains our data: First, using the mean-field expression for x, in 
order to reproduce the experimental dip, the peak in \ at T c has to be at least an order of 
magnitude larger than that suggested by magnetization measurements 1 . Second, from our 
resistivity data we obtained H c2 (T), and hence we calculated £,o(T), l(T), an d lastly, X' L . 
Our values of X' L also are unable to explain the magnitude of the drop of A below T/v — a 
conclusion also shared by Gammel et al. from their SANS data- 

Since the CN model does not take into account the effect of the SVP on the superfluid 
density, we neglect it when we convert to superfluid density p s : First we assume AA(T) 
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FIG. 1: (O) Temperature dependence of the penetration depth AA(T) from ~0.1 K to 13K. Inset: 
AA(T) below 6.5 K. The arrows show features at 6 K (AFM), 2.3 K (WFM) and 0.45 K (SVP). 

follows a power-law temperature dependence at low temperatures from the combination of 
gap-minima observed in non-magnetic borocarbides and the increased pairbreaking as Er 
spins disorder. Consequently, we set \i ow (T) = A(0)(1 + bT 2 ) with b = 0.036 K -2 from 
Ref. 0. Next we offset Xi ow until it matches the data at 1.3 K, the local minimum in AA in 
the vicinity of 1.5 K. Finally we convert AA to p s in Fig. [21 

The superfluid data leads to some important observations. First, the data in the param- 
agnetic phase (T > T/v) fit the theoretical curve based on an isotropic superconducting gap 
(solid squares), and not that based on nodes. The solid line shows a superfluid calculation 
based on a d x 2_ y 2 order parameter. Second, the superconductivity is only slightly depressed 
in the AF phase below T N . The best fit to data above T N (solid squares) is obtained when 
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FIG. 2: (O) Experimental superfluid density p 8 (T) = [A 2 (0)/A 2 (T)] calculated from AA(T) data 
in Fig. ^ from 0.1 K to T c . Solid squares = Calculated p s (T) assuming paramagnetic phase from 
T=0 to T c . Note that 1/tAo=24, 26 or 28 give virtually the same theoretical curve. Solid line: 
Theoretical curve for d x 2_ y 2 order parameter. The arrow denotes Neel temperature at 6 K. 

A(0) = 740 A — here we assume paramagnetism from T— to T c , neglecting AF order, with 
parameter 1/tAq = 24. The paramagnetic curve is almost unchanged if one uses l/VAo=26 
or 28. We see that the paramagnetic curve fits the data above T/v, and overestimates the 
data in the AF phase. We will show below that, as one crosses T/v from above, the AF 
phase leads to increased or decreased superfluid density depending on the combined effects 
of the following three factors: (1) the AF molecular field, which decreases the magnitude of 
the superconducting gap and hence decreases the superfluid density, (2) freezing out of spin 
fluctuations, leading to decreased pairbreaking and hence increased superfluid density, and 
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FIG. 3: Chi-Nagi model calculation for the superconducting gap A. Solid line = BCS temperature- 
dependence. (O) Paramagnetic gap. Solid squares = Incorporating AF phase below Tjv, for 
different values of 1/tAq. Note that when 1/tAo=28, superconductivity is fully recovered. 

(3) scattering from non-magnetic impurities, which reduces the suppression of the gap by 
the molecular field. The parameters of ErN^P^C are such that these three effects result in 
a slight decrease in superfmid density below T N . 

Fig. El shows the calculated superconducting gap amplitude A(T), in the presence of 
Hq(T), for various values of l/rA , as well as the paramagnetic curve. The normalized 
paramagnetic gap (open circles) agrees excellently with the BCS gap (solid line). Tunneling 
measurements^ also show that A(T) follows the BCS curve above T^. Next, as shown in 
Fig. in the AF phase, as l/rA increases, superconductivity is gradually recovered, as 
evidenced by the increase of A(T) (solid squares). Fig. 0] shows p s (T) for various values of 
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FIG. 4: Superfluid density p s {T)= [A 2 (0)/A 2 (T)] from 0.1 K to T c . Solid circles = data. AF-phase 
fitting curves for l/rA =24 (□), 26 (A), 28 (x). 

1/tAq. Notice that when l/rA =26, (1) p s decreases only slightly at T N , and (2) A is only 
slightly depressed below the BCS value, in agreement with tunneling dat a 31 ' 32 . This value of 
1 / t Ao corresponds to a mean free path (mfp) of 45 A. Though this curve still overestimates 
the experimental data below Tjy, it at least fits the data better than the paramagnetic curve. 

To see if this value of 1/tAq is reasonable, we take A(0) = 1.83T C for ErNi 2 B 2 G2&, for 
which the BCS coherence length ^ BCS = ^f_p/7rA(0) = 470 A, where vp = 3.6 x 10 5 m/s 
is taken from band-structure calculations^ 3 - for LuNi 2 B 2 C and YNi 2 B 2 C Using the relation 

# c2 (0)= 0.693T c (^)| r , (16) 
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and dH C 2/ (1T\t c {H || c) = —2.67 kOe/KA, we obtain the coherence length 



Finally, using the relation^ 

^ = 0.85(£ A SC5 1/2 , (18) 

we obtain the mfp 1=42 A. On the other hand, from the resistivity value just above T c , 
p(T c *)=5.8 /if2-cm, we get 1=56 A. These two values agree well with the value of 45 A 
calculated from l/VAo=26 obtained earlier, implying that this particular value of the non- 
magnetic scattering rate needed to explain our p s data is consistent with other measurements. 
Note that this value of mfp calculated from H c2 data does not depend on the exact value of 
T c . Also, the prefactors 0.693 and 0.85 in the above relations are for materials in the dirty 
limit. Here I < £q, so our ErNi 2 B 2 C sample may be considered as "quasidirty" . It is puzzling 
that our sample has a high T c and be considered quasi-dirty, yet this is consistent with the 
results of other papers. Also, in this sample the non-magnetic scattering rate (1/t) is at 
least two orders of magnitude larger than the effective magnetic scattering rate (l/r 2 ), 
thus the mfp value is largely determined by 1/t. Our mfp value, however, is smaller than 
the 90 A obtained from resistivity measurements just above T c in Ref. |35j The CN model is 
thus able to explain our superfluid density data, both qualitatively and quantitatively. Our 
data, in agreement with others, also shows that AF order coexists with superconductivity 
below T N . 

It is also interesting to note that according to Ref. 0, a near-exact cancellation of spin- 
fluctuation and molecular-field effects occur at a critical value of N(0) J c ^~1.0 x 10~ 3 , where 
J c f =1 | gj — 1 | is the conduction-electron local (/) spin exchange. For the case of ErNi 2 B 2 C, 
we obtain N(0)J c f = 1.7 x 10~ 3 , which explains the small change in pairbreaking at T^. 

We turn next to an alternative explanation for the change of p s at Tjy. Ramakrishnan 
and VarmaF^ predicted that for materials with a nested FS, since the peak in susceptibility 
and the joint density of states (defined as the difference between the susceptibility in the 
superconducting state and the normal state) occur at the same Q-value, one should expect 
an increase in pairbreaking at Tjy. Conversely, a non-nested FS will give rise to decreased 
pairbreaking at T/y. Two-dimensional angular correlation of electron-positron annihilation 
radiation measurements show that only one out of the three FS sheets in LuNi 2 B 2 C possesses 
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FIG. 5: p s (T) from 0.1 K to T c . (O) Data. (A) Calculated AF-phase curve for l/rA =26, (□) 
Calculated AF-phase curve incorporating WFM. Note that Twfm here is chosen to be 3 K. 

nesting properties, thereby accounting for the propensity for magnetic ordering found in the 
other magnetic members of the RE nickel borocarbides 27 . Also, Dugdale et ali^ estimated 
that the fraction of the FS that would be able to participate in nesting is only 4.4%. Contrast 
this with CN model, which assumes perfect 1-D nesting. Hence the increased pairbreaking 
due to partial nesting on one FS sheet is partially compensated by decreased pairbreaking 
by the other two sheets, resulting in only a slight increase in pairbreaking at T^. 

As temperature further reduces below T/v, the theoretical curve in Fig. El overestimates 
the experimental curve below 3 K (~0.3T C ). This is due to additional pairbreaking effect of 
the ferromagnetic moments in the WFM phase, which shows up as a small peak near 2.3 K 
(see Fig.P). The small dip in superfluid density shows that this WFM slightly depresses, but 
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does not completely destroy, superconductivity, demonstrating the coexistence of WFM and 
superconductivity. We model this WFM by including a temperature-dependent magnetic 
impurity scattering rate 1/t^ fm = l/20{ gj - 1) J (J + 1)(1 - T/T WFM ) U ' (with the same 
value of gj and J as before), and adding this to the previous effective magnetic scattering 
rate, i.e. l/r 2 totei = l/r 2 e// + 1/t^ fm when T < T WFM . The pre-factor 1/20 arises from 
the fact that one out of every 20 spins contributes to the WFM 3 ^, giving rise to a weak 
magnetization. The temperature dependence {1 — T /T WFM ) U ' is analogous to the molecular 
field formulation. We obtain v' « 2 from Jensen's calculation 37 or Choi and Canfield's 
data 3 ^. Note that p s already starts to flatten out at 3 K, consistent with neutron-scattering 
data, which shows that this weak ferromagnetic component already shows up at 3 K— . 
Hence we choose Tw FM =3 K in this WFM calculation. Fig. |5] shows p s when one accounts 
for WFM. The calculated p s does flatten out below 3 K, but does not increase below 2.3 K 
as the data did. There is as yet no direct measurement of superconducting gap amplitude 
at this temperature range, though our model predicts a drop in A there. One may need 
to include the effect of inelastic spin-fluctuation scattering in this low temperature region, 
which was ignored by the CN model. 

In conclusion, we present in-plane penetration depth data of single- crystal ErNi 2 B 2 C 
down to 0.1 K. The small increase in pairbreaking at T/v can be attributed to the interplay 
between the effects of the AF molecular field and spin-fluctuation scattering. It could also 
be due to the combined effects of non-perfect nesting on one piece of the FS and non-nesting 
on other pieces of the FS. The increased pairbreaking at Twfm is modelled by a magnetic 
impurity scattering parameter, and both magnetic orders coexist with superconductivity. 
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